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We analyze the sigma meson mass and width together with the pion-pion scattering parameters 
in the context of non-local chiral quark models with wave-function renormalization (WFR). We 
consider both non-local interactions based on the frequently used exponential form factor, and on 
fits to the quark mass and renormalization functions obtained in lattice calculations. In the case 
of the sigma properties we obtain results which are less dependent on the parameterization than in 
the standard local NJL model, and which are in reasonable agreement with the recently reported 
empirical values. We also show that the inclusion of the WFR tend to improve the description of 
^ , the 7T-7T scattering parameters, with the lattice inspired parameterization providing the best overall 

y-^ ■ results. Finally, we analyze the connection of the non-local quark models discussed here with Chiral 

' Perturbation Theory, and present the model predictions for the low energy constants relevant for 

\Q ' n-TY scattering to 0(4) in the chiral expansion. 



I. INTRODUCTION 



Although much effort has been made in trying to predict low energy hadron observables directly from QCD, one 
is still far from reaching this goal due to the extremely complex non-perturbative behavior of the theory in that 
regime. In such a situation it proves convenient to turn to the study of effective models. For two light flavors it is 
believed that QCD supports an approximate SU(2) chiral symmetry which is dynamically broken at low energies, 
and pions play the role of the corresponding Goldstone bosons. A simple scheme including these properties is the 
well known Nambu— Jona-Lasinio (NJL) model [3], proposed more than four decades ago. The NJL model has been 
widely used as an schematic effective theory for QCD @, d, 3, allowing e.g. the description of light mesons as 
fermion-antifermion composite states. In the NJL model quarks interact through a local, chiral invariant four-fermion 
^} coupling. Because of the local nature of this interaction, the corresponding Schwinger-Dyson and Bethe-Salpeter 
equations become relatively simplified. However, the main drawbacks of the model are direct consequences of this 



locality: loop integrals are divergent (and therefore have to be regulated somehow), and the model is nonconfining. As 
a way to improve upon the NJL model, extensions which include nonlocal interactions have been proposed (see Ref. Q 
. £h ' an d references therein). In fact, nonlocality arises naturally in quantum field theory and-particularly, in several well 
established approaches to low energy quark dynamics, as ejp the instanton liquid model jo] and the Schwinger-Dyson 
resummation techniques Q. Lattice QCD calculations @, [UBS also indicate that quark interactions should act over 
5r ! a certain range in the momentum space. Moreover, it has been argued that nonlocal extensions of the NJL model do 
not show some of the above mentioned inconveniences of the local theory. Indeed, nonlocal interactions regularize the 
model in such a way that anomalies are preserved [ll| and charges are properly quantized, the effective interaction is 
finite to all orders in the loop expansion and therefore there is not need to introduce extra cutoffs [12] , soft regulators 
such as Gaussian functions lead to small next-to-leading order corrections [l3j], etc. 

In the present work we will reconsider non-local models adopting as the basic ingredient a reliable description of 
the quark propagator as given from fundamental studies, such as lattice QCD. In this sense, it should be noticed that, 
except for Ref. flB] . most of the calculation performed so far using non-local chiral quark models have neglected the 
wave function renormalization in the propagator (See e.g. Refs. ( 1 6l. Tl7l. Hii lisl l20lp . Recent lattice QCD calculations 
suggest, however, that such renomalization can be of the order of 30 % (or even more) at zero momentumfl, 0, licj ]. 
Moreover, these calculations also show that the quark masses tend to their asymptotic values in a rather soft way. 
Thus, it is of importance to perform a detailed study on the incorporation of these features in this type of models, 



'Electronic address: Santiago.Noguera@uv.es 
tElectronic address: scoccola@tandar.cnea. gov. ar 



2 



and analyze their role in the prediction for different hadronic observables. The lagrangian we will use is the minimal 
extension which allows to incorporate the full momentum dependence of the quark propagator, through its mass and 
wave function renormalization. Using this lagrangian we explore which are the implications for some pion and sigma 
meson properties originated by changes in the quark propagator. In particular, we present here results for the sigma 
meson mass and width, and for the pion-pion scattering parameters. Studying these scattering parameters close to 
the chiral limit we are also able to obtain predictions for some of the low energy constants of the Chiral Perturbation 
Theory (xPT) Lagrangian 

The present article is organized as follows. In Sec. II we present the model lagrangian and the formalism necessary 
to derive some selected pion and sigma meson properties. In Sec. Ill we discuss different ways to obtain the model 
parameters and compare the resulting quark propagators with available lattice data. In Sec. IV we present and 
discuss the predictions of the model for the selected parametrizations, paying special attention to the role played 
by the incorporation of the wavefunction renormalization and by the difference in the quark interaction momentum 
dependence. In Sec. V we analyze the connection of the non-local quark models described here with xPT, and present 
the predictions for the corresponding low energy constants relevant for 7r-7r scattering to 0(4) in the chiral expansion. 
Finally, in Sec. VI our main conclusions are summarized. 



II. THE MODEL 
A. Effective action 

Let us begin by stating the Euclidean action for the nonlocal chiral quark model in the case of two light flavors, 

Se = J d A x t^(x)(-ip + m c )i)(x)-^Y j a (x)j a (x)-jp(x)jp(x) j . (1) 

Here m c is the current quark mass, which is assumed to be equal for u and d quarks. The nonlocal currents j a (x),jp(x) 
are given by 



Ja(x) = I d 4 z g{z) ip (x + T a ip (x - |) . 



r — r 

Here, T a = (1, rysf) and u(x ) d v(x) = u(x')d x v(x)—d x >u(x')v(x). The functions g(z) and f(z) in Eq.(|2]), are nonlocal 
covariant form factors characterizing the corresponding interactions. The four standard quark currents, j Q (x), require 
the same g(z) form factor to guarantee chiral invariance. The new term, j>(x)jp(x), is self-invariant under chiral 
transformations. The scalar-isoscalar component of the j a (x) current will generate the momentum dependent quark 
mass in the quark propagator, while the "momentum" current, jp(x), will be responsible for a momentum dependent 
wave function renormalization of this propagator. For convenience, we take the same coupling parameter, Gs, for 
the standard chiral quark interaction and for the new jp(x)jp(x) term. Note, however, that the relative strength 
between both interaction terms will be controlled by the mass parameter >c p introduced in Eq.((2]). We have choosen 
the relative sign between these terms in order to have a real value for x p for the case in which the wave function 
renormalization Z (p) (explicitly defined in Eq. ([10] below) is less than 1. In what follows it is convenient to Fourier 
transform g(z) and f{z) into momentum space. Note that Lorentz invariance implies that the Fourier transforms g(p) 
and f(p) can only be functions of p 2 . 

In order to deal with meson degrees of freedom, one can perform a standard bosonization of the theory. This is 
done by considering the corresponding partition function Z = JVi/jVip exp[— Se], and introducing auxiliary fields 
<7i(x), 02(2;), 7?(x)), where 01,2 (x) and 7?(x) are scalar and pseudoscalar mesons, respectively. Integrating out the quark 
fields we get 

Z = j Vox Va 2 Vtt ey^[S E ° s ] , (3) 



where 



S h E os = - In dctA+ -i- / -jpL [ ai (p) cri (_p) + ^( p ) . *(_ p ) +<72 ( p) a2 (-p)] . (4) 



2GsJ (2tt) 
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The operator A reads, in momentum space, 

A(p, p')={-j> + m c ) (2tt) 4 (p - p') + .g [ti (p' - p) + ^ ■ *(p' - p)\ 

+ f 



p+p'\ i>+f 

2 X n 



0"2 ip' ~P), 



(5) 



At this stage we assume that the oi l2 fields have nontrivial translational invariant mean field values 012, while the 
mean field values of the pseudoscalar fields 7Tj are zero. Thus we write 



0-2(2;) = >c p 02 + 5o 2 (x) 
tt(x) = 5tt(x) 

Replacing in the bosonized effective action, and expanding in powers of the meson fluctuations, we get 

obos oMFA 1 oq ua d , 



(6) 
(7) 
(8) 



Here the mean field action per unit volume reads 



cMFA 
°E 



-2N r 



d A p 



trln^ 1 ^)] 



2G.c 



2 -2 

2G S 



where the quark propagator in the mean field approximation T>o{p) is given by 

Zip) 



-p + M (p) 



with 



The quadratic terms can be written as 



Z(p) = (1 - o 2 fip))- 1 
M(p) = Zip) (m c + a x gip)) 



^quad 1 

E ~ 2 



d 4 P 



[G a ( P 2 ) 5a(p) 5a(-p) + G a ,( P 2 ) 5a' (p) Sa'(-p) + G^p 2 ) 5tt( p ) ■ 5x(-p)] 



where the a and a' fields are related to 01 and 02 by 

5a = cos 5a\ — sin 6 5a 2 
5a' = sin 9 8a\ + cos 8 5a 2 



(9) 



(10) 



(11) 



(12) 



(13) 
(14) 



and the mixing angle 9 is defined in such a way that there is no a — a' mixing at the level of the quadratic action. 
The function G^ip 2 ) introduced in Eq. (fl2"|) is given by 



G 7r ip 2 ) = — -8N i 



9 2 (q) 



Z{q + )Z(q-) 
D(q+)D(q-) 



q- +M(q+)M(q-)] 



with q = q ±p/2 and D(q) = q 2 + M 2 (q), while for the a — a' system we have 



where 

G ai ai ip 2 ) 

G<j 2 <j 2 ip 2 ) 

G ai a-2 ip 2 ) 



G 



1 

1 

gS + 

8N C 



ip 2 ) = 



G ai a x (P 2 ) + G a2<T2 ip 2 ) 



T \ [G^aAP 2 



G<7 1 a 1 ip 2 ) — Ga^oiip 2 ) 



n 2 



(15) 



(16) 



8iV c 



d 4 q 



9 2 il) 



Z(q+)Z(q-) 



d 4 q 



(2tt)4 * ™D(q+)D{q 

d 4 q 2.2. Zjq+)Z(q-) 
(2tt) 4 q 1 {q) Diq+)D(q-) 
Z(q+)Z{q-) 



J- [q+ .q- -M{q+)Miq-)} 



gii)fiq) 



Diq+)D(q 



(q + ■q-)-M(q + )Miq-) + 
^-q-[q-Miq + ) + q+M(q-)] 



(g+) 2 (g-) 2 - (g+ -g-f 
2q 2 



(17) 
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B. Mean field approximation and chiral condensates 

In order to find the mean field values <J\,2, one has to minimize the action S^ FA . A straightforward exercise leads 
to the coupled gap equations 

<j\ — 8N C Gs J ^2^yr d(p) ~D{p) = 

/j4 2 Z(v) 

wr^ p fip) W) =0 - (18) 

Now the chiral condensates are given by the vacuum expectation values (qq) = (uu) = (dd). They can be easily 
obtained by performing the variation of Z MFA = cxp[— S E IFA ] with respect to the corresponding current quark 
masses. This expression turns out to be divergent. Thus, as customary, we regularize it by subtracting its value for 
non-interacting quarks. We obtain 



• qq} 4iVc i (2tt)4 V D(p) pZ + ml) 



C. Meson masses and quark-meson coupling constants 

The meson masses can be obtained by solving the equation 

G M {-m 2 M ) = . (19) 
In the case of the a — a' system the mixing angles is given by 9{—m 2 a where 

tm2 ^>= oJgr£Iw (2o » 

Finally, the on-shell meson-quark coupling constants gu q q are given by 

, _ dGM{p) 



9Mqq = ^Mqq(~ m M) 



dp 2 



(21) 



Note that due to the mixing, in the scalar meson channel the corresponding vertex has two components. Thus for 
aqq vertex we have 



o i i + if 

Vcrqq = 9 aqq 1 + 9aqq~^ 



(22) 

p 



where 



gl qq= 9aqq C0S9 ] 9^ q \ = 9aqg SUX9 (23) 

D. Pion weak decay constant 

By definition the pion weak decay constant f„ is given by the matrix element of the axial current (x) between 
the vacuum and the renormalized one-pion state at the pion pole: 

(0\A;(0)\« b (p)}=iS ab P »U . (24) 

In order to obtain an explicit expression for the axial current, we have to "gauge" the effective action Se by 
introducing a set of axial gauge fields A't (x) . For a local theory this "gauging" procedure is usually done by performing 
the replacement 

d^d^+ l - lb T-A^{x) . (25) 
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In the present case — owing to the nonlocality of the involved fields — one has to perform additional replacements in 
the interaction terms. Namely, 



ip(x-z/2) -> W A (x,x - z/2) ip(x-z/2) 
^{x + z/2) -> ipi(x + z/2) W A (x + z/2,x) 



(26) 



Here x and z are the variables appearing in the definitions of the nonlocal currents (see Eq.j2])), and the function 
Wa(x, y) is defined by 



W A (x,y) = P exp 



(27) 



where s runs over an arbitrary path connecting x with y. 

Once the gauged effective action is built, it is easy to get the axial current as the derivative of this action with 
respect to At (x), evaluated at A^x) = 0. Performing the derivative of the resulting expressions with respect to 
the renormalized meson fields, we can finally identify the corresponding meson weak decay constants. After a rather 
lengthy calculation we obtain 



five Qitqq 



with 



F (p 2 



8N C 



(27T) 



D(q+)D(q-) 



q~ + M(q+)M(q-)] 



(28) 



(29) 



It is important to notice that the integration over the path variable s appearing in this calculation turns out to be 
trivial and, thus, the result path- independent. In the chiral limit the expression Eq. ((28|) has a rather simple form 
given by 



M (0) 



(30) 



which connects with the Goldberger-Treiman relation. 



E. The decay width of the Sigma meson 

To obtain the decay amplitude of the a meson into two pion we need to calculate 

(2tt) 4 5\q + qi + q 2 ) S ab G a ^(q 2 , q 2 ,q 2 ) 



SS% os 



5a{q)5Tr a (qi)6Tr b (q 2 ) 



(31) 



where the meson fields are assumed to be already renormalized. In terms of the unrenormalized fields and taking into 
account the <j\ — a 2 mixing we have 



where 



Ga^(q 2 ,ql,ql) =G<j q q(q 2 ) G vqq (qf) G nqq (q 2 ) G anv (q 2 , qf, q 2 ) 



G alT7r {q ,q 1 ,q 2 ) = G ai7rn (q ,q 1 ,q 2 ) cos% ) - G CT27r7r (g ,q 1 ,q 2 ) sin% 



(32) 



(33) 



and the expressions of the unrenormalized u\ and a 2 coupling constants to two 7r can be obtained by expanding T to 
third order in the fluctuations. We get 



G ai ^{q 2 ,q 2 ,q 2 2 ) = -IQN C J 



f h+k 2 \ f k + kA f k + k 2 \ Z(k)Z(k x )Z{k 2 ) 
9 \ 2 J 9 V 2 ) 9 { 2 J D(k)D(h)D(k 2 ) X 

M(k) ki ■ k 2 - M(fcj) k ■ k 2 - M(k 2 ) k-fa- M(k) M(fci) M{k 2 ) 



(34) 



G a2 ^{q 2 , qj,q 2 ) = - 



8iV c 



/ fci + fc 2 \ / fc + fci \ ( k + k 2 ^ z(fc)z(fci)z(fc 2 ) x 



(27T) 



D{k)D{k 1 )D{k 2 ) 



x kjk-h + kj k ■ k 2 + (fci + k 2 ) ■ [M(fci) M(k) k 2 - M{k 2 ) M(h) k + M{k 2 ) M(k) k x 



(35) 



where k\ = k + qi and k 2 = k — q 2 and q 2 = (q\ + q 2 ) 2 . Similarly for a' we have 

G a '^{q 2 ,ql,ql) = G a > q q(q 2 ) G^qqiqf) G^qqiq 2 ) G^^iq 2 , qf, qf) (36) 

where 

G al7 , 1T {q 2 ,ql,q 2 )^G (Jl ^{q 2 ,ql,ql) sm6(q 2 ) + G a2 ^(q 2 , qf, q 2 ) cos% 2 ) (37) 
In terms of gu-n-n = GMTrir('m 2 VI , m 2 , rn^ ) the M = a, a' width reads 



v 3 9m™ L 4m2 

2 167rmM V 

F. it-it scattering 

In general, the total amplitude for the ir-ir scattering process can be expressed as 

A(iT a (q 1 ) + 7T l3 (q 2 ) -> 7r 7 (g 3 ) + it 6 (q 4 )) = 5 a p5 l5 A(s,t,u) + 5 ai 5p 5 A(t, s,u) + 6 aS S^A(u,t, s) (39) 

where 

s = (<Zl+<72) 2 ; t = { qi -q s ) 2 ■ u = {q 1 -q i ) 2 (40) 

Within the present model, this amplitude gets two contributions. One corresponds to the box diagram and the other 
to the scalar meson pole diagram. Thus, 

A(s,t,u) = A box (s,t,u) -gi qq G m™( s > m l' m l) G m( s ) ( 41 ) 

where 

A box (s,t,u) =gi q q [J(s,t,u) + J(s,u,t) - J(u,t,s)] (42) 

and 

J {s,t,u) = i [J box (q 1} q 2 ,q 3 ) + J box -q 3 , -92)] (43) 

with 

d i k fk + kA (k + k 2 \ fh + k 13 \ fk 2 + k 13 \ Z(ki)Z(k)Z(k 2 )Z(k 13 ) 



Jfcoz (?i, 92, 93) = 16iV c / — g I — - — I .9 I — ^ — I 



2 7 D{ki)D{k)D{k 2 )D{k l3 ) 
{ [fci • fc + M(fci)M(fc)] [fc 2 • fcis + M(fci)M(fci 3 )] - [fci • fc 2 + M(fci)Af(fea)] [fc • fc 13 + M (fc)Af(fc 13 )] 

+ [fcj ■ fc 13 + M(*i)M(fc 13 )] [* • fc 2 + Af(*)M(fe)] } (44) 

where k x = k + qi , k 2 = k - q 2 , fc i3 = k + qi - q 3 . 

It is customary to define the scattering amplitudes of defined isospin 

T° = 3A(s, t, u) + A(t, s, u) + A(u, t, s) 
T 1 = A(t,s,u) - A(u,t,s) 

T 2 = A(t 7 s, u) + A(u, t, s) (45) 

In terms of these amplitudes the scattering lengths a\ and slope parameters b\ are defined by the partial wave 
expansion at low q 2 

J' dx P t (x) T z (s, t, u) = q 21 [a\ + b\ q 2 + ...) (46) 
where Pe(x) is the Lagrange polynomial of order I. 
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III. DETERMINATION OF THE MODEL PARAMETERS 



In this section we present in some detail the procedure used to determine the model parameters as well as the form 
factors g(q) and f(q) which characterize the non-local interactions. 

In our first model (scenario SI) we use exponential functions to model the non-local interactions. These are well 
behaved functions which have been often used in the literature (see e.g. [I^.[l7l.fT§.[l9|) to define g(q). Here, we also 
use such form for f(q). Thus, for SI we have 

g(p)=cM-P 2 /K) ; f(p) = cxp(-p 2 /A?) (47) 

Note that the range (in momentum space) of the nonlocality in each channel is determined for the parameters Ao and 
Ai, respectively. From Eq. (jXTJ we obtain 

M (p) = Z (p) [m c + cti exp( V/Aq)] 

Z( P )= [l-cr 2 exp(V/A?)] _1 (48) 

We fix the values of m c and < qq > 1 ^ 3 to reasonable values m c = 5.7 MeV and < qq > 1 / 3 = —240 MeV determining 
the rest of the parameters so as to reproduce the empirical values f v = 92.4 MeV and = 139 MeV, and Z(0) = 0.7 
which is within the range of values suggested by recent lattice calculations [5 fioj] . 

For the second parametrization we follow Ref.(l4l|. where a parametrization based on a fit to the mass and renor- 
malization functions obtained in a Landau gauge lattice calculation was used. Such parametrization is 

M (p) =m c + a m f m {p) , 

Z{p) = l + a z f z (p) , (49) 



with 



Up) = i + (pVAlf 2 ; Up) = [i + ( P 2 /Aj)] ^ , (so) 



-5/2 



where the analytical form of f m (p) has been proposed in Ref.0|. The analytical form of f z (p) is chosen in order 
to guarantee the convergence of the integrals. Some alternative parametrization of this type suggested from vector 
meson dominance of the pion form factor can be found in Ref.|22l|. In terms of the functions f m (p) and f z (p), and 
the constants m c , a m , a z the form factors g(q) and f(q) are given by 



g(p) 



l + a z a m f m (p) - m c a z f z (p) 
l + ctzfz (p) OL m - m c a z 



1 + OtzJzKP) 



and the mean field values are 



n fTlc^-z 

1 + a, 



°2 = t— — 52 

l + a z 

The parameters for this second model (scenario S2) are determine as follows. As before we take Z(0) = 0.7 and fix 
Ao and Ai in such a way that the functions f m (p) and Z (p) agree reasonable well with lattice results of Ref.[j|. Next 
we fix m c and a m in order to reproduce the physical values of m w and / w . The resulting parameters are m c = 2.37 
MeV, a m = 309 MeV, and with A = 850 MeV and Ai = 1400 MeV. 

Finally, in order to compare with previous studies where the wavefunction renormalization of the quark propagator 
has been ignored we consider a third model (scenario S3). In such scenario we take Z(p) = 1 and exponential 
parametrization for g(p). Such model corresponds to the "Model II" discussed in Ref.[l^|. from where we take the 
parameters corresponding to < qq > 1 / 3 = —240 MeV. 

The values of the model parameters for each of the chosen scenarios are summarized in Table [H In FigQ] we 
compare the quark mass function f m (p) and renormalization function Z(p) as obtained from our three scenarios with 
data extracted from the lattice results of Ref . @ ■ The main reason for comparing f m {p) (instead of M (p)) is that 
analyzing lattice data from different groups using Landau gauge fixing[3, [l3|, and also results for M (p) obtained 
by each group using different inputs, we observed that the resulting functions f m (j>) are very similar in spite of the 
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differently looking M(p). On the other hand, the renomalization functions Z (p) are much less sensitive to the choice 
of lattice parameters, and in fact the two lattice groups [1, provide similar results. We observe that the functions 
fm (p) and Z (p) for scenario SI, based on exponential functions, decrease faster than the lattice data. For scenario 
S2, however, they go to zero as (p 2 ) -3 / 2 and (p 2 ) -5 / 2 , respectively, following the lattice data in a closer manner. 
Finally, in the case of S3 the exponential decrease of f m (p) is even faster than that of Si . 

IV. NUMERICAL RESULTS 

In this section we present and discuss our numerical results. In Table U we give the results for the mean-field 
properties, together with the pion and sigma masses and decay parameters. As it can be seen in this table, while 
for the exponential parameterizations (i.e. SI and S3) the empirical values of f n and ra^ are consistent with a 
quark condensate which lies within the range of the usually quoted phenomenological values — {qq) 1 / 3 — 200 - 
260 MeV (23[ l2~ij the scenario S2 leads to a value of the chiral condensate somewhat above such range. On the 
other hand, the corresponding current quark mass is quite smaller than those obtained for the scenarios SI and S3. 
This issue deserves some comment. The chiral condensate, as well as the current quark masses, are scale dependent 
objects. In particular, the phenomenological values quoted above for the condensate correspond to a choice of the 
renormalization scale p = 1 GeV. In the parametrization S2 some parameters have been determined so as to obtain 
a good approximation to the lattice mass renormalization function Z(p), a quantity which also depends on the 
renormalization point. In particular, we use the function Z(p) obtained in Ref.Q where the renormalization scale 
has been chosen to be p = 3 GeV. One might wonder whether the fact that this renormalization point differs from 
the one usually used to quote the values of the condensate can account for the fact that the S2 prediction is outside 
the empirical range. If one assumes that this difference is also responsible for the rather low value of m c this can be 
investigated in the following way. To leading order in the chiral expansion the current quark mass and the condensate 
are related by the Gell-Mann-Oakes-Renner (GMOR) relation 

fl ml = 2 < qq > m (53) 

where rh = (m u + m<i) /2 The validity of GMOR to that order is well justified by the low energy behavior of the irir 
scattering amplitudes [25[. Using that, according to Ref . (26| . rh runs from 5.5 MeV at the scale p = 1 GeV to 4.1 MeV 
at p = 2 GeV we expect that a typical value of (<?<?) 1/3 = -240 MeV at p = 1 GeV will run to (qq) 1/3 = -270 MeV 
at p = 2 GeV [13] ■ Lattice calculations provide an independent determination of quark masses and qq condensate 

is mm 

m£f (2 GeV) = 4.3 ± 0.4 stat +1A MeV (54) 

U.4 sys 

{qq) (2 GeV) = - (265 ± 5 stat ± 22 sys MeV) 3 (55) 

which confirms the p = 2 GeV values given above. Note that since these two lattice calculations are not connected, the 
quoted values imply a verification of the GMOR relation. Since the GMOR relation is well satisfied by our lagrangian 
model [lij], and in all our scenarios f„ and are fitted to their empirical values, it is clear that the quality of the 
description of the quark condensate and the current quark mass are closely related. Thus, a further running up to 
p — 3 GeV implies that the current quark mass must be scaled by a factor of the m(2 GeV)/m(3 GeV) = 1.11. This 
value is rather different from the factor 1.81 obtained from the ratio between the lattice result at p = 2 GeV and the 
value of m c for the scenario S2 given in Table [H This clearly indicates that possible ambiguities related to the choice 
of renormalization point cannot fully account for the rather high value of the condensate for the scenario S2. In fact, 
using the above mentioned factors to reescale the value — {qq) 1 / 3 — 326 MeV quoted in Table U down to p = 1 GeV 
we get —{qq) 1 / 3 — 284 MeV which is about 10% above the empirical upper limit. A possible way to reduce the value 
of the quark condensate in S2 is to reduce the parameter Ao. For Ao ~ 600 MeV we can obtain values for the quark 
condensate and quark masses which are within the phenomenological bounds. 

The mass and width of the sigma meson display some dependence on the parametrization. However, such dependence 
is smaller than the one found in the local NJL model[3l|. The obtained values for the masses are somewhat larger than 
the recently extracted empirical values 478^23 ± ^7 MeV Q and 390t3gMeU [13| while the widths are compatible 
with the experimentally reported values 324^ 2 ± 21MeU Hi and 282±^MeU (33|. 

The situation concerning the a' meson deserves some comment. In general, for the non-local models under con- 
sideration the quark propagators develop a series of poles in the complex plane. In Euclidean space, such poles can 
be purely imaginary (as in the NJL model which only has one pole of this type) or fully complex. The existence of 
these poles implies the appearance of "pinch points" [IB] in the calculation of the meson two-point functions. The 
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external momentum for which the first of such "pinch points" appears is given by p pp = 25, where Si is the imaginary 
part of the first pole of the quark propagator. From this point on the functions G in Eq.(17) do in general develop 
an imaginary component related to the unphysical decay into qq pairs, which is usually associated with the lack of 
confinement. In some cases, depending on the regulator and/or parametrization, one can find a prescription for the 
integration path along the complex plane such that this imaginary component cancels outfit. [l7l . l34| . It is clear, 
however, that the corresponding results turn out to be prescription dependent and, unless the meson pole appears 
no far above p pp , not very reliable. For this reason, in this work we take the point of view that p pp marks the limit 
of validity of our model. For the three scenarios under consideration we have found p pp to be about 1 GeV, which 
appears to be a reasonable scale for a low energy effective model of QCD. As for the a' channel we have verified that 
no pole corresponding to a meson of this type appears below that scale. 

We turn now to the low-energy parameters for 7r — tt scattering. These parameters have been matter of much 
attention in the recent past years. In particular, recent results on K14 decays (35l [36| have led to an improved 
phenomenological determination [3?], EH] of the threshold parameters for S-, P-, D- and F- waves. Our results for the S 
and P waves are displayed in Table HU while those corresponding to D and F waves in Table Mil Since the calculation 
of sigma pole contributions include off-shell quantities it is not possible to perform a clear and unique separation 
between a and a' contributions. Thus, only the sum of such contributions is given. In general, reasonable estimates 
indicate that a' contributions represent only a few percent of this total value. The phenomenological values extracted 
in Ref.[38| are also indicated. In comparing our results with these values one should keep in mind that the present 
model does not incorporate pion loops, and hence there is still room for improvement. Finally, for comparison, in 
Tables ITT1 and IIIII the existing results for the local SU{2] NJL model [39l.l40| are given. Results obtained in alternative 
QCD-based quark models can be found, e.g. in Ref.|4ll| 

We analyze first the results corresponding to the S- and P-waves. Let us recall that to leading order in the chiral 
expansion the corresponding length and slope parameters can be obtained from the Weinberg amplitude 

A(a,t,u) = ^fi (56) 

which leads to the predictions 

^ a° ■ m„ = -4 a 2 ■ m v = b° ■ ml = -2 b 2 ■ m\ = 6 a\ ■ ml = (57) 

Since our three different scenarios lead to the same values of f„ and m^ the predictions for these five scattering 
parameters are expected to be quite similar. In fact, results in Table HU confirm this, although those of S2 are in 
slightly better agreement with empirical data. This is particularly interesting in the case of a 2 ,, which results from 
a rather strong cancellation between box and sigma contributions. In order to be more sensitive in the comparison 
between scenarios, we also give in Table ITT1 the combination of the S-wave isospin and 2 parameters 2a§ + 7o,g which 
vanishes in the chiral limit. We observe that in all scenarios the correction goes in the right direction. Moreover, in 
the case of S2 its magnitude is larger providing therefore a better description of the experimental result. Another 
way to improve on the discrimination between the different parametrizations of our model is to consider corrections 
up to q e order in the expansion Eq. (|46|) . Thus, we calculate the parameters cf and d\ corresponding to the q 4 and 
q e corrections, respectively. We observe that in each partial wave the exponential interaction produces scattering 
parameters which decrease rather fast with the power of q 2 . On the other hand, the scenario S2 predicts coefficients 
which are of the same order of magnitude in each partial wave. 

We consider now the scattering lengths and slope parameters for D- and F-waves displayed in Table IIIII These 
results, together with the scattering lengths and slope parameters of S- and P-waves given in Table HTl complete all 
cases for which there are phenomenological determinations available. For SI and S3, we observe that the signs of the 
parameters are correctly predicted, except for 6° m S3. The absolute values for the scattering lengths are off by a 
factor between 1.5 and 2.5, whereas the slope parameters fail by one order of magnitude. On the other hand, the 
scenario S2 gives the right sign and order of magnitude in all cases, deviating only by a factor 3 in the worse case, 63. 

From the previous results we can conclude that although the exponential interaction might be able to reproduce 
the scattering lengths parameters rather well the description of higher power coefficients is, in general, expected to 
be less accurate as the power in q 2 increases. This is particularly so for the higher partial waves. On the other 
hand, the momentum dependence of the scenario S2 seems to be better adapted for the description of the higher 
power parameters. In fact, the only case where this scenario gives a worse result than the exponential ones is in the 
prediction for b\, where a strong cancellation between the box and sigma contribution takes place. 

Comparing scenarios SI and S3 we can observe the effect on the scattering parameters of taking into account the 
wave function renormalization. Except for the parameters listed in Eq. lj57|l . we observe that as the power in q 2 
increases the associated parameters obtained in scenario S3 decrease faster than in scenario SI. We can conclude that 



10 



the effect of the wave function renormalization term goes in the right direction, even if this effect is less important 
than the one produced by the difference in the momentum dependence of the interactions. It should be noticed that 
our scenario S3 is very similar to the model used in Ref.j42j. In fact, in both cases the wave function renormalization 
is not included, exponential parameterizations are used and the values of and f„ are fitted. The difference comes 
from the way in which the third parameter of the model is determined. In Ref.[42j] the rather sensitive value of a| 
was used, while here we choose to fix the chiral condensate. 

In our scenarios which include wavefunction renormalization we have fixed Z(Q) = 0.7. As it can be seen in FigHJ 
however, for small values of p the errors in the lattice results are rather large. Thus, Z(0) is not well constrained 
by lattice calculations. In order to test the sensitivity of our results to this kind of uncertainties we have reduced it 
to Z(Q) = 0.6, and considered the scenario S2 for two alternative situations. In the first case we varied the model 
parameters so that /„■ and remain at their empirical values, while in the second case we kept the model parameters 
fixed. In both cases we found that most of our results change by less that 10%, the most notorious exception being 
the 7T7T scattering length a| which changes about 15 %. It is interesting to note that in the second case the pion mass 
and decay constant, as well as the chiral condensate, get reduced. We obtain = 138.7 MeV, f n = 91.2 MeV and 
— < qq > 1/3 = 323 MeV. 



V. COMPARISON WITH CHIRAL PERTURBATION THEORY 



In the previous section we have focused our attention on the ability of our quark model to reproduce the phe- 
nomenological tt-tt scattering parameters. An alternative point of view (see, for example Refs. [39L lia. H. I4H I4H l47l ) 
is to consider the quark models as the generators of the pion Chiral Perturbation Theory (xPT) Lagrangian [21 1. 
%PT describes the low energy physics of pions in a universal way, once the order in the momentum and chiral break- 
ing expansion (i.e. the order in the chiral expansion) is specified. Different scenarios for quark models will lead to 
XPT Lagrangians with different values of the so-called low energy constants (LECs) . In this section we analyze the 
connection between our quark scenarios and the xPT Lagrangian up to the fourth order in the chiral expansion. 

To perform this connection we introduce the pionic Lagrangian 



where 



£ = £2+ U , (58) 



£2 = ^(d^U* d»U + tf X + X ] U) , (59) 
£4 = £1 (dptf d»U) 3 + £2 (d^ d v U) d v U) 2 + £ 3 ( X U) 2 + £4 (0„x d^U) 2 + ... , (60) 



and 



f/ = exp(V^J ; X = m 2 (; i ) . (61) 

Note that among all possible 0(4) terms only those relevant for tt-tt scattering to that order have been explicitly 
given. To the order we are working, the parameters / and m can be related with the predicted values for f n and 
through 

/ - / ( 1+ (x) >< ')- (62) 

n4=m a (l + 2feV/ 3 ) . (63) 



Using (|62| and l(63|) . we can express the scattering parameters resulting from Eq.(|58[i in terms of the £4 coupling 
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constants and the m v and / ff values as follows 

= 32tt I y r ' 7 V /* 



'l + if^l [40 4+40 ^ + 104 + 14 4] 




m„ a 9 = 



■5s(t) {2A + 4) (64) 

As already mentioned, the Lagrangian lj58|) is valid up to fourth order in the chiral expansion, therefore it can be fully 
equivalent to our quark model scenarios only when they are treated to the same order of approximation. Thus, to 
extract the LECs defining the pionic Lagrangian from the values of the scattering parameters, f„ and obtained 
in each of our quark scenarios we should analyze the values of these parameters for small values of m c . In fact, we 
have verified that close to the chiral limit the scattering parameters display, as a function of (m T // T ) 2 , the quadratic 
behavior expected from Eqs.(|64|). From the determination of the corresponding linear and quadratic coefficients it 
is possible to obtain the numerical values of LECs 4 It should be noticed that this procedure for obtaining ii is 
completely equivalent to the bosonization of the quark Lagrangian followed by a covariant gradient expansion (see 
Ref.[39l| for the application of such method to the NJL model). At this stage we are connecting our quark model at 
the one loop level to the pionic Lagrangian at the tree level. Our next step is to make the connection with the xPT 
Lagrangian. The main difference between the Lagrangian (|58| and the xPT Lagrangian is that our ii parameters are 
finite and no pion loop contribution is present. The scattering parameters in xPT [21| include pion loop contributions, 
and are written in terms of renormalized LECs i\. As expected, Eqs. (|64|) coincide with the ones obtained from xPT 
if the corresponding pion loop contributions are neglected. In this approximation the coupling constants ii can be 
identified with the l\ constants at some given renormalization scale [i. 

Eqs. f64|) imply several relations between the scattering parameters. We focus on two of them 



Test! = m T (2 a° - 5 a 2 ) + ml (~a\ -b° + h£) = J m j 

\ Z ' \ 16ir 4 /* 



using JM 



^ using X PT (65) 



5 2 3 n f using JJ1 . . 

Test2 = - m^ a 2 + —m^b x - m n aJ, = i l 7tt ■ pT (66) 
2 10 L 16tt 4 /* 450 using x^J- 

Obviously these two relations vanish when we use our pionic Lagrangian (|58)l at the tree level. Therefore, a non- 
vanishing value obtained for these two quantities in any other calculation must be originated by loop corrections or 
by higher order terms in the chiral expansion. As indicated in Eqs. (|65|66[) . in the case of the xPT Lagrangian both 
relations have corrections from pionic loops. On the other hand, the deviation from zero of Testl and Test2 when 
evaluated using the scattering parameters obtained in our quark scenarios at the physical value of m^ is originated by 
higher order terms in the chiral expansion. In TablelTVlwe show the results for these two relations in our scenarios. Also 
indicated are the xPT Lagrangian results, which corresponds to the pion loop contribution of the order (m T /27r/ T ) . 
From this table we observe that the quark scenarios previously studied give results for Testl and Test2 which are of 
the same order of magnitude than the pion loop contributions. This implies that the studied quark models include 
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higher order contributions (i.e. 0(6) or higher) which are as important as the chiral loops. The effect of these higher 
order contributions is more important for the scenario S2, due to its different behavior for large momenta. 

In Table V we give the li values corresponding to our different scenarios. It is interesting to note that in the case of 
t\ the listed values result, in all cases, from an important cancellation between the box and the sigma contributions. 
For I2 only box contribution is present since no scalar meson contribution is possible [48j]. Also given in Table V 
are the values of the renormalized LECs l\ (/i) obtained[3~3] in the framework of xPT at some particular values of 
renormalization scale ^[Ilj. We observe that the sign and order of magnitude of the most accurately known LECs 
l\ and l\ are well reproduced for small values of [i. In fact, in the case of the scenarios SI and S2 the agreement 
is remarkable good for /1 around 2 which is a reasonable scale since we have integrated out degrees of freedom 
from below the sigma mass. In the case of the LECs l\ and 1%, even though it is not so good in the case of S2, the 
agreement is still acceptable given the existing uncertainties in the determination of the empirical values. Finally, as 
a reference, some typical values obtained within the local NJL model taken from Ref.[39] are also listed in Table V. 

VI. CONCLUSIONS 

In this work we have analyzed the sigma meson mass and width together with the ir-ir scattering parameters in 
the context of non-local chiral quark models with wave-function renormalization (WFR) term. We have considered 
two types of momentum dependence for the quark interactions. The first one (scenario SI) is based on the frequently 
used exponential form factors. The second one (scenario S2) corresponds to a fit to the mass and renormalization 
functions obtained in lattice calculations Q, and gives rise to a softer momentum dependence (e.g. at large momentum, 
the quark mass decreases as (p 2 ) -3 / 2 instead of exponentially). In order to test the influence of the WFR, we also 
considered a third scenario S3 which corresponds to an exponential interaction but where this renormalization is 
absent. 

Our results for the sigma mass are relatively stable, ranging from 552 MeV for S2 to 683 MeV for S3. We observe 
that the coupling between the scalar term of the standard chiral interaction and the new scalar term associated with 
the WFR term reduces the value of the lower sigma mass, as it must be expected. Comparing the S3 and SI results 
we observe a reduction of a 10%, while in the case of the S2 interaction there is a further reduction of 10% originated 
by the softer momentum dependence of the interaction. The width of the sigma follows the same reduction as its 
mass, as one goes from one scenario to another. These results are less dependent on the parameterization than in the 
standard NJL model. The predicted mass and width are reasonable close to the recently reported empirical values 

turn. 

Regarding the ir — ir scattering parameters, we have compared our results with the phenomenological determination 
made in Ref.jH]. Although the existence of the chiral limit relations, Eqs. lf57|) . for the S- and P-wave scattering length 
and slope parameters reduces the sensitivity of these parameters to the choice of the different quark interactions, we 
have been able to discriminate between these interactions by going to higher order in the momentum expansion or 
to higher partial waves. We conclude that although the exponential interaction is able to reproduce the scattering 
lengths parameters rather well the description of higher power coefficients turns out to be, in general, less accurate 
as the power in q 2 increases. This can be clearly seen in the case of higher partial waves. On the other hand, the 
momentum dependence of the scenario S2 seems to be better adapted for the description of the existing empirical 
data. Comparing the predictions of the scenarios based on exponential interactions, SI and S3, we observe that the 
presence of the WFR term tends to improve the results, even though its effect is less noticeable that the one produced 
by the difference on the momentum dependence of the interactions. 

Finally, we have analyzed the relation of our quark scenarios with the chiral Lagrangian up to 0(4) in the chiral 
expansion. In particular, we have obtained predictions for the low energy constants li involved in it-it scattering 
within our scenarios. The procedure we followed, using the scattering parameters, is equivalent to the standard 
method of bosonization followed by covariant gradient expansion. Our predicted values for li are in relative good 
agreement with the values for the renormalized l\ constants defined in the xPT calculations 12 ill for a /i value about 
2m T . They are also in the range of values obtained in the NJL model calculation of Ref. [23]. We have been able 
to define combinations of the scattering parameters which allow to discriminate between higher chiral corrections 
(0(6) or higher) and pion loop corrections. We observe that the higher order corrections included in our non-local 
quark model calculations at physical are of the same order that the pion loop corrections not considered in this 
work. The effect of such corrections in our predictions for the mesonic observables is an issue that deserves further 
investigation. 
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Table I: Model parameters and results for some alternative parameterizations. 







SI 


S2 


S3 




MeV 


5.70 


2.37 


5.78 


G s Aq 




32.03 


20.82 


20.65 


A 


MeV 


814.42 


850.00 


752.2 


yep 


GeV 


4.18 


6.03 




Ai 


MeV 


1034.5 


1400 




0\ 


MeV 


529 


442 


424 


<T2 




-0.43 


-0.43 




M(0) 


MeV 


375 


311 


430 


Z(0) 




0.7 


0.7 


1.0 


-<OT> 1/3 


MeV 


240 


326 


240 



MeV 139 139 139 

g^qq 5.74 4.74 4.62 
jw MeV 92.4 92.4 92.4 



Ilia 


MeV 


622 


552 


683 


(0) 
Qcrqq 




5.97 


4.60 


5.08 


Q W - 

ijoqq 




-0.77 


-0.26 






MeV 


263 


182 


347 
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Table II: 7T — 7T scattering parameters for S and P waves. 





contribution 


SI 


S2 


S3 


NJL 
Ref.[39] Ref.140] 


Empirical 
Ref.[38] 


f \ 

(m w ) x 


box 


-1.536 


-1.279 


-1.618 








a 


1.718 


1.470 


1.798 








Total 


0.182 


0.191 


0.180 


0.18 0.19 


0.223 ± 0.009 


(m 3 ) x 6g 


box 


0.114 


0.117 


0.114 








a 


0.116 


0.146 


0.107 








Total 


0.230 


0.263 


0.221 


0.22 0.27 


0.290 ± 0.006 


(m„.j x c 


box 


— U.UUoO 


U.Uzoo 


— U.UU I D 








a 


U.U41Z 


U.uooo 


U.UoUz 








total 


U.Uozo 


n AQ07 


U.Uzzo 






(m n ) x a 


box 


u.uuuo 


u.uoo 


U.UUU4 








(7 


U.UUo I 


n mo 
u.uiy 


U.UUoi 








Total 


o nono 

u.uuyz 


(1 nor 
U.UOO 


O OAK £ 

U.UUOO 






/ \ 2 
(m-rr) x a 


box 


-0.6851 


-0.5790 


-0.7170 








a 


0.6404 


0.5346 


0.6721 








Total 


-0.0447 


-0.0444 


-0.0449 


-0.046 -0.044 


-0.0444 ± 0.0045 


(m 3 ) x 6g 


box 


-0.053 


-0.049 


-0.051 








a 


-0.031 


-0.034 


-0.033 








Total 


—0.084 


—0.083 


—0.084 


—0.091 —0.079 


—0.081 ± 0.003 


fro 5 \ y 


box 


u.uuou 


n DD7S 
u.uu / o 


U.UUOZ 








(7 


0.0042 


0.0056 


0.0034 








Total 


0.0121 


0.0134 


0.0116 






( m 7 ) X dn 


box 


-0.0005 


-0.0006 


-0.0005 








(7 


-0.0006 


-0.0011 


-0.0004 








Total 


-0.0011 


-0.0017 


-0.0009 






77171- X (^(Zo ~r ' ^0 J 




U.UDz 


U.U / z 


U.U40 


U.U4 U.U/z 


n i 4- n aq^ 


(m% 10 ) x aj 


box 


25.1 


23.9 


24.7 








a 


10.5 


11.3 


11.1 








Total 


35.7 


35.2 


35.7 


37 34 


38.1 ± 0.9 


(m* 10 3 ) x b\ 


box 


5.56 


4.60 


5.34 








a 


-2.10 


-2.85 


-1.72 








Total 


3.45 


1.75 


3.62 




5.13 ±0.15 


(ml 10 3 ) x c\ 


box 


0.21 


-2.70 


0.15 








a 


0.38 


0.63 


0.25 








Total 


0.59 


-2.06 


0.40 
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Table III: Scattering lengths and slope parameters for D and F waves. 







contribution 


C 1 

SI 


CO 

Si 


C 9 

So 


NJL 
Rpf T^Ql Rpf T401 


Empirical 

SXxiL.l OOJ 


(ml 


10 4 ) x ai 


box 

(7 


9.71 

a on 
4.ZU 


9.76 

0.0 / 


9.93 

'iAA 
0.44 










lotal 


lo.y 1 


1 K A Q 


1 Q Q7 

16.6 1 


iq 7 i a 7 
lo. / 10. / 


lo. 66 ± U.oO 


(ml 


10 4 ) x b° 2 


box 

(7 


0.98 

1 Ofi 
— l.Zo 


0.85 

o on 
— z.zu 


1.04 
— U.oO 










Tn+ n 1 


u.ou 




U.IO 






(ml 


10 4 ) x a\ 


box 
o 


-2.74 
a on 

4.ZU 


-2.95 
o.o / 


-2.43 

'iAA 
0.44 










lotal 


1.40 


O 70 
Z. / Z 


1.U1 


11 Q O 
1.1 O.Z 


/ifi 4- n oc; 
Z.40 ± U.ZO 


(ml 


10 4 ) x 6l 


box 

(7 


0.08 

1 Ofi 
— l.Zo 


0.07 

o on 
—z.zu 


0.13 
—U.oO 










TVif ^1 

lotai 


1 1 Q 
— i. iy 


9 1/1 
— Z. 14 


n 73 

— U. f 6 




— o.oy it u. io 


(ml 


10 5 ) x a\ 


box 
a 


0.82 
1.82 


1.15 
3.09 


0.7 
1.2 










Total 


2.65 


4.24 


1.9 




6.05 ±0.29 


(ml 


10 5 ) x b\ 


box 
a 


0.06 
-0.70 


0.0 
-1.6 


0.07 
-0.40 










Total 


-0.64 


-1.6 


-0.33 




-4.41 ±0.36 



Table IV: Results for Testl and Test2 defined in Eqs. (|65|66[) as obtained in our quark scenarios (S1,S2 and S3) and Chiral 
Perturbation Theory to 0(4) (xPT). The results obtained using the empirical values of Ref.[38] (Empirical) arc also given. 





SI 


S2 


S3 


X PT 


Empirical [38] 


Testl x 10 2 


-1.2 


-2.5 


-1.1 


1.5 


0.40 ± 4.4 


Test2 x (m 4 10 4 ) 


0.1 


-3.3 


0.006 


1.6 


3.21 ± 1.4 



Table V: Values of £i obtained in our different scenarios. The xPT values of £\ as a funtion of n are obtained from Ref.[33|. The 
last two columns corresponds to the NJL predictions from Ref. for two different constituent quark mass: M — 220,264 MeV. 

Non Local Quark Model X PT (^(/-P) NJL 



SI S2 S3 fi = m p i i = 2m- w i i = m 1! M — 220 M = 264 

£1 x 10 3 -1.39 0.26 -2.07 -4.0 ± 0.6 -1.9 ±0.6 -0.4 ± 0.6 -0.63 -2.28 

£2 x 10 3 6.46 6.41 6.51 1.9 ±0.2 6.2 ± 0.2 9.1 ±0.2 6.29 6.18 

£3 x 10 3 -2.3 -4.1 -1.1 1.5 ±4.0 -1.8 ±4.0 -4.0 ±4.0 -8.50 -3.48 

£4 x 10 3 17.2 20.3 15.0 6.2 ±1.3 19.1 ± 1.3 27.9 ± 1.3 22.73 12.16 
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Figure 1: (Color online) fm(j>) (Upper panel) and Z(p) (lower panel) for various parametrization as compared with Lattice 
results of Ref. @] 



